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Abstract 

In this paper the thermodynamics of thin ferromagnetic films are stud- 
Q"\ ied in the framework of the classical Heisenberg model with uniaxial 

anisotropy and long-range dipole interaction. The dependence of the 



m 

00 

o 
m 



order of the reorientation transition in these systems on the number 
of layers and on the distribution of anisotropies is investigated using 
both mean field theory and Monte-Carlo simulations. 
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Thin magnetic films with a ground state magnetization perpendicular 
to the surface can undergo a phase transition into a phase with in-plane 
magnetization at a temperature Tr < Tq, the Curie temperature of the film, 
due to the competition of dipole interaction and surface anisotropy [1]. This 
reorientation transition has already been investigated theoretically in the 
framework of the two-dimensional anisotropic Heisenberg model with dipole 
interaction using classical [2] and quantum mechanical [3] mean field theory, 
renormalization group theory [4], spin wave theory [5], and Monte-Carlo 
simulations [6, 7]. While all calculations for a monolayer except [6] show a 
first order reorientation transition, the calculations in [3] show a second order 
transition for a system with L = 3 layers and vanishing anisotropy in the 
middle layer. In a recent paper on the order of the reorientation transition 
Moschel et al. [8] proposed that the order depends on the distribution of 
anisotropics through the layer. 

In this paper we present ground state aspects and detailed phase dia- 
grams obtained from both mean field theory and Monte-Carlo simulations 
in a classical anisotropic Heisenberg model with long range dipole interac- 
tions. The system consists of L two-dimensional layers and it is described 
by the Hamiltonian 

+ "J J2 r rj 3s * ' S 3 ~ 3r ij 5 ( S * ' T ij)( T ij ■ S J')> 

ij 

where the Sj = (sf,s^,sf) are classical Heisenberg spins at position = 
(rf,ry,rf) on a simple cubic lattice and r^- = — rj. J is the nearest- 
neighbour exchange coupling constant, K\ is the uniaxial anisotropy in layer 
A = 1 . . . L, and D = n^g 2 '47ra 3 is the strength of the long range dipole 
interaction on a lattice with lattice constant a. We consider ferromagnetic 
coupled systems with J = 1 in this paper. 

One open question concerning this model is the occurence of stripe do- 
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mains in the ground state of the monolayer (L = 1). Taylor et al. [2] showed 
that for certain values of the dipole interaction D and large anisotropy K\ 
there are stripe domains with all spins pointing into the indirection, but 
they were unable to extend their calculations into the range where D is 
small. Furthermore they proposed that for small D the two-dimensional 
system should be in a ferromagnetic state, ie. the stripe domain width N^ s 
should diverge for a finite value of D. This uncertain result called for further 
investigation. 

We analytically calculated the ground state stripe domain width Nq S 
and the ground state energy Eq S , with the result 

Nq S = N e^D~ + O ( e ~^) 

and 

Egs = ~Ki -2J + d(^-- + O [e~^ , 

with the constants iVo = 0.871026 and the two-dimensional dipole sum [9] 
0o = 9.0336. Note that Nq S is finite for all finite values of D, although it 
grows exponentially fast when J/2D gets large. 

In the limit of large Nq S the ground state properties of the system can 
be determined exactly: We find that a system with only one layer (L = 1) 
always has a first order phase transition with a free energy independent 
of the orientation of the spins when K\ = 6.77522Z?. When L = 2, the 
transition is first order only when the anisotropics in both layers are equal 
(K\ = K2 = 6.52962.D), while it is second order when K\ and K2 are 
different. For L = 3 the situation is more difficult, as we only get a first 
order transition when K\ = A'3 = 6.52920.D and K2 = 6.28402Z?. These 
ground state results indicate that the occurence of a first order transition is 
a rather special case in this model. 

We investigated the finite temperature regime with mean field theory 
and found finite regions in parameter space where the system undergoes a 
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first order reorientation transition. In figure 1 the phase diagram for L = 2 
layers, J = 1, D = 0.017 and K\ + Ki = 0.25 is depicted. At the phase 
boundaries the in-plane-component (xy) and the ^-component (z) of the 
total magnetization become unstable. The lines of second order transitions 
are dashed and the lines of first order transitions are solid. The arrows shall 
symbolize the orientation of the spins. 

In the Monte-Carlo simulations we considered systems with N x N x 
L spins and periodic boundary conditions in the in-plane directions. The 
precise method of the simulation is described in [7]. The results of the 
Monte-Carlo simulations are in very good agreement with the mean field 
calculations: For L = 1 we only find first order reorientation transitions, 
while for L = 2 the transition is of second order when K\ differs enough from 
K2- To illustrate these results, in figure 2 the histogram of the polar angle 
$ of the total-magnetization vector in layer 1, is depicted. When the 

transition is of second order (o), the distribution has one maximum around 
while at a first order transition (o) the distribution has two maxima, 
one at t?i = and a second at i?i = tt/2 due to hysteresis effects. The 
differing results in [6], where a second order transition is found for L = 1, 
are very likely due to short simulation times and an insufficient evaluation 
method. While in [6] a typical run involved 4 x 10 5 MC sweeps, we performed 
2 x 10 6 sweeps and additionally used a global spin flip, which reduced the 
autocorrelation time by a factor of approx. 200. A detailed account of our 
work will be published elsewhere. 

Acknowledgement: One of the authors (A. H.) would like to thank 
A. Moschel for fruitful discussions. This work was supported by the Deutsche 
Forschungsgemeinschaft through Sonderforschungsbereich 166. 
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Figure captions 



Figure 1: Mean field phase diagram for L = 2 layers, J = 1, D = 0.017 and 
K\ + K2 = 0.25. At the phase boundaries the in-plane-component (xy) and 
the z-component (z) of the total magnetization become unstable. The lines 
of second order transitions are dashed and the lines of first order transitions 
are solid. The arrows shall symbolize the orientation of the spins. 

Figure 2: Histogram h(&\) of the polar angle of the total magnetization 
vector in layer 1, as obtained from MC simulations, for 32 x 32 x 2 spins. 
J = 1, D = 0.017, and K x + K 2 = 0.25, with K L = 0.05(c) and K x = 
0.125(o). The lines are guides to the eyes. 
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Figure 2 
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